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ABSTRACT 

I propose an analysis method, based on spin-spherical harmonics and spherical Bessel 
functions, for large-scale weak lensing surveys which have source distance information 
through photometric redshifts. I show that the distance information can significantly 
reduce statistical errors on cosmological parameters; in particular, 3D lensing analysis 
offers excellent prospects for constraining the equation of state of the vacuum energy 
which dominates the energy density of the Universe. I show that the ratio of pressure 
to energy density could be determined to an accuracy of ~ 1% or better. Having 
distance information also offers significant advantages in the control of systematic 
effects such as the intrinsic alignment of galaxies. The case for obtaining photometric 
redshifts is therefore compelling. A signal-to-noise eigenmode analysis of the modes 
shows that the modes with highest signal-to-noise correspond quite closely to ignoring 
the redshift information, but there is significant extra information from a few radial 
modes. These modes are generally long-wavelength, suggesting that useful information 
can be gleaned even if the photometric redshifts are relatively inaccurate. 

Key words: cosmology: observations - gravitational lensing - large scale structure, 
galaxies: formation 



1 INTRODUCTION 

As a direct probe of the mass distribution, gravitational lens- 
ing is an excellent tool for cosmological parameter estima- 
tion, complementing microwave background studies. Many 
of the cosmological parameters have been determined well 
from the microwave background, most recently and most ac- 
curately by the WMAP satellite (Spergel et al. 2003). There 
are some parameters which are difficult to constrain from 
the microwave background alone, the most significant ex- 
ample being the equation of state parameter {w — p/p(?) 
of the vacuum energy. An accurate determination of w and 
its time-dependence is a major goal in the post- WMAP era, 
as it can put constraints on the nature of the vacuum en- 
ergy. The direct mass dependence also offers advantages over 
studies of large-scale structure, which have to make assump- 
tions about, or deduce, the relationship between galaxy and 
mass distributions on large scales (Verde L. et al. 2002). One 
of the most useful manifestations of gravitational lensing by 
intervening matter is the alignment of nearby images on the 
sky. Detection of dark matter on large scales through such 
'cosmic shear' measurements has recently been shown to be 
feasible, and cosmological constraints from the first detec- 
tions have now been made (e.g. Bacon, Refregier & Ellis 
2000; Van Waerbeke et al. 2000; Brown et al. 2002a). To 



date, the cosmic shear analyses have concentrated on using 
essentially only two-dimensional shear information; in sur- 
veys where (photometric) redshifts of the sources are known, 
they are used to determine the redshift distribution only, and 
are not used individually. Recently, though, 3D information 
has begun to be employed, mostly through what I would de- 
scribe as 2iD analysis: the sources are divided into slices at 
different distances, based on their photometric redshifts, and 
a 2D analysis is done on each slice. This can be valuable in 
detecting the presence of clusters of galaxies (Wittman et al. 
2001; Wittman et al. 2002; Taylor et al. 2003). At a statisti- 
cal level, Hu (1999) has shown that there is some extra infor- 
mation on cosmological parameters which can be gained by 
dividing the sample this way. The extra information depends 
on the parameter under consideration: for the amplitude of 
the matter power spectrum, there are significant gains when 
the source population is split into two, but little is gained by 
further, fine subdivisions. For w, the gains are much larger, 
as the constraints from 2D lensing are weak. Finally, Taylor 
(2001) has shown how the lensing equations can be inverted 
to deduce the gravitational potential directly, and this offers 
exciting possibilities of reconstruction of the 3D mass den- 
sity field from shear and convergence data (Bacon & Taylor 
2002). Hu & Keeton (2002) have shown how this inversion 
can be done by matrix methods for discrete shear data in 
cells. 
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The purpose of this paper is to demonstrate how the full 
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3D information about the shear can be used. In particular, 
I investigate how source photometric redshift information 
can be used in a statistical analysis to determine the mass 
power spectrum and the equation of state of the vacuum 
energy. In particular, wc arc intorostod in how much extra 
information on cosmological parameters can be gleaned from 
a 3D shear map, to see whether the reduction in statistical 
errors justifies the extra observing time required to obtain 
photometric redshifts for the sources. 

Reducing statistical errors is, of course, not the only 
reason for wishing to have photometric redshifts; weak sheax 
studies are beset by a number of systematic errors, and these 

may eventually dominate the error budget as surveys im- 
prove in scope and quality. The most obvious of these is 
the redshift distribution of the sources - ignorance of the 
precise distribution is already a source of significant error 
in current surveys (see e.g. Brown et al. 2002b). At a dif- 
ferent level, the analysis of cosmic shear usually makes the 
assumption that the source galaxy ellipticities are uncorre- 
lated, and this may not be true in detail, due to tidal effects 
during the galaxy formation process (Heavens, Refregier & 
Heymans 2000; Croft & Metzler 2000; Crittenden et al. 2001; 
Catelan, Kamionkowski & Blandford 2001; Mackey, White 
& Kamionkowski 2002; Jing 2002). The level of this effect is 
not well known, but it is thought to be large enough to be 
an important source of systematic error in low-rcdsliift Icns- 
ing surveys (Brown et al. 2002b; Heymans & Heavens 2003). 
By its nature, tidal effects are important for galaxies which 
are close together in three dimensions, so photometric red- 
shifts can be used to remove or downweight nearby galaxy 
pairs in shear correlation estimates, for example Heymans 
& Heavens (2003), King & Schneider (2002). These stud- 
ies show that the systematic effects of intrinsic alignments 
can be essentially completely removed, at the expense of an 
increase in shot noise. Other physical effects which could af- 
fect weak Icnsing studies can also be removed this way, for 
example the effects of source clustering if one uses number 
counts as a measure of magnification (Broadhurst, Taylor 
& Peacock 1995). Furthermore, with many (typically 5) im- 
ages of galaxies available from a photometric redshift lensing 
survey, there is scope for better shape measurement, inde- 
pendent lensing studies and so on. 

We see that there is a powerful case for obtaining photo- 
metric redshifts for lousing surveys on the basis of removing 
systematic errors. Having obtained distance information for 
the sources, it seems sensible to use it in the statistical analy- 
sis, as it must reduce the error bars on the parameters which 
one wants to estimate. In this paper, I propose a method, 
based on spherical Bessel functions and spin-spherical har- 
monics, for analysing such surveys. 

I find that the 3D information can reduce the error on 
the amplitude of the matter power spectrum by a factor 
~ 2, which broadly speaking compensates for the extra time 
taken to obtain photometric redshifts. The error on the vac- 
uum equation of state is encouraging: with 3D information w 
could be determined to an accuracy of < 1%. These results 
assume that other cosmological parameters are known accu- 
rately, most plausibly from microwave background studies. 

The paper is arranged as follows. Section 2 describes 
the method, section 3 discusses the expected errors, section 
4 performs a signal-to-noise analysis of the modes, which 



turns out to be very informative, and section 5 presents the 
conclusions. 



2 METHOD 

2.1 Spherical coordinates and spherical Bessel 
functions 

Spherical coordinates are natural for the description of fields 
on the full sky, but there are reasons for using them even 
if the sky coverage is relatively small. The survey specifi- 
cation is normally a combination of sky coverage (specified 
in terms of angular coordinates 6,(fi), and depth (related 
to r). In addition, photometric errors introduce errors in r, 
and the lensing potential is related to the Newtonian po- 
tential by a radial integral. Having said this, the spherical 
harmonic expansion becomes cumbersome at high £, and at 
some point a Fourier description on the flat sky becomes an 
attractive approximation to make. We choose to expand ra- 
dially in spherical Bessel functions, as these (when combined 
with spherical harmonics) are eigenfunctions of the Laplace 
operator, which leads to a very simple relationship between 
the coefficients of the gravitational potential and those of 
the overdensity field, since these are related by Poisson's 
equation. 



2.2 Transformation of scalar and shear fields 

The spherical harmonic transform of a scalar field /(r) is 

defined by 



4/ 



rf=^r/(r)j,(fer)yr*(") 



(1) 



where ji{z) is a spherical Bessel function, Y a spherical 
harmonic, fe is a wavenumber, ^ is a positive integer and 
m = —£, . . .£. The shear field, comprising 71 and 72, can be 

written as components of a tensor, and expanded in terms 
of tensor spherical harmonics, as in Taylor (2001), or spin- 
spherical harmonics sYl^ as in Okamoto & Hu (2002). We 
will use the latter notation, generalised to 3D: 



71 (r) ± 172 (r) = 



El 



{£-2)\ 



(2) 



dkk ■yim{k)±2Yf"{n)je{kr). 



In principle the coefficients -ytmik) could depend on whether 
the -|- or — sign is taken, but they are the same, and are 
related to the transform of the lensing potential 4>{r) by 



7<!' 



.ik) = 2 



1 (£ + 2)1 



i£-2)\ 



(jitmik) 



(3) 



(Taylor 2001). The lensing potential is related to the gravi- 
tational potential $(r) by a radial integral (e.g. Bartelmann 
& Schneider 2001) 



(r) = 



dr' 



fK{r)-fK{r') 



fK[r)fK{r') 



<I.(r'). 



(4) 



fK{r)dij) is the dimensionless transverse comoving sep- 
aration for points separated by an angle dip. The 
Robertson- Walker metric may be written ds^ = c^dt^ — 
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R^{t)\dr^ + f^{T■)d■lp'^^^, and fxij-) takes the values 
sinr, r, sinhr for curvature values fc = 1, 0, — 1. 

3>(r) is related to the overdensity field 5(r) = 5p{v)/p 
by Poisson's equation 



2a{t) 



(5) 



where Q,m is the matter density parameter, Hq is the Hub- 
ble constant and o(t) = 1/(1 + z) is the scale factor. 5 itself 
is not a homogeneous field, because it evolves with time, 
and hence with distance from the observer through the light 
travel time. In the linear regime, we can describe the growth 
by a universal growth rate g{t) (see, for example Abaza- 
jian & Dodelson (2002) for quintessence models), so that 
5 = S^g{t), where (5°(r) is the overdensity field, extrapo- 
lated using linear theory to the present day, and which is 
homogeneous. Thus #(r) = ^°g{r)/a{r). Note that for an 
Einstein-de Sitter universe, g{r)/a{r) = 1, and that for the 
'concordance' model of matter and vacuum density param- 
eters Om = 0.3, SIl. = 0.7, g{r)/a{r) declines from 1.28 at 
high redshift to unity today. 

The transform of the lensing potential is therefore re- 
lated to that of the present-day potential field by 



1 f°° 

(t>em{k) = - dfc'fc'V(fc,/c')*L(fc' 
Jo 

= - / drfK 
^ Jo 



(6) 



where 

Ve{k,k' 



f 

Jo 



dr'je{k'r 



{r)jt{kr) 

fK{r) - fK{r') 



9{r') 
a{r') ' 



(7) 



These equations allow us to relate the coefficients of the 

shear field, 7^m(fc) to the underlying linear, present-day over- 
density field, Sim{k), since Poisson's equation implies 

*L(fe) = -^^^5L(fc). (8) 

2.3 Discrete estimator 

Tiie data which we have to hand are estimates of the shear 
field at 3D positions in space, so it makes some sense to 
transform the data directly. The radial coordinate is gen- 
erally not known accurately, but is typically given by a 
photometric redshift which has an error which may be 
ffz ~ 0.02 — 0.1 or more. We denote this radial coordinate 
by s, and the true coordinate by r. They are related by a 
conditional probability, which we will take to be a gaussian 



p{s\r)ds = 



27ra"; 



exp 



{Zs 



2al 



(9) 



where Zr,s are the redshifts associated with distance coordi- 
nates r and s. Note that Oz can vary with z. 

The estimate of the shear is obtained from the complex 
ellipticity (Bartelmann & Schneider 2001) 



e = -I- 27 



(10) 



where e" is the intrinsic source ellipticity, whose mean is as- 
sumed to be zero, so 7 = e/2 is an unbiased, albeit noisy, 
estimate of the shear field. The shear estimates have covari- 
ance properties 



(7a7b> = (7a76> + 7(6066*) 



(11) 



for galaxies a and 6. We denote the variance of ehy a^, typ- 
ically < 0.1 (Hudson et al. 1998). The second term isnon- 

zero if galaxies are intrinsically aligned (see e.g. Heavens, 
Refregier & Heymans 2000). I ignore this term, as with pho- 
tometric information, its effect can be removed by modifying 
the analysis to remove close pairs in three dimensions (Hey- 
mans & Heavens 2003; King & Schneider 2002). Thus we 
have noisy estimates 7 of the shear field at imprecise posi- 
tions with distances s given by photometric redshifts. 

An obvious set of quantities to use for a harmonic de- 
scription of the data is simply to do a discrete transform 
of the measurements which are to hand. We consider the 
following quantities: 



±ge, 



galaxies g 

±gem{k) are estimates of the quantities 
±gemik) = / d^sn{s)f±{r)je{ks) ±2Y£"'(n), 



(13) 



where n(r) is the number density, and 7± = 71 ± 172. It is 
straightforward to include a non-uniform weighting of the 
galaxies if desired. This may improve the errors on parame- 
ter estimation, but is not explored in this paper. 

Note that in this equation the shear field is evaluated at 
the true position r, whereas the estimate (12) involves the 
distance s estimated from the photometric redshift. These 
are related by the conditional probability p{s\r), leading to 
an average value of the expansion coefficients 



±gem{k) = 



d^s J drp{s\r) 
n{s)j±{r)ji{ks) ±2Y«'"(n). 



(14) 

With the photometric redshift smoothing, n(s) is heavily 
smoothed, so we can approximate it by the smoothed num- 
ber density at r, ris (r) . For deep surveys, the angular cluster- 
ing is small, so we can ignore clustering, and approximate 
this by the average number density, which normally sepa- 
rates into a radial part and an angular selection: 



ns{r) = n{r)M{n) 



(15) 



Normally M = (unobserved sky) or M = 1 (in survey), 
although more complicated forms are possible, if some parts 
of the survey are partially sampled. 

If we denote the spin-spherical harmonic transform of 
ns(r)7±(r) by ±hem{k), i.e. 

ns(r)7±(r) = y'^^ J dkk'^ ±hem{k)je{kr)±2Yl"{n) (16) 

we find that (14) may be written 
±gem{k) = ^ J dsf](:{s)dndrp{s\r) 

'ie'm'ik)j'e{k'r)jeiks) 



/ dk'k'^ ±he 



±2Y,T'{n) ±2yr*(ft)M'(n). 
It is convenient to define angular mixing matrices by 



(17) 
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iW^r = I dn ±2F,7 (n) ±2Fr*(fi)M2(n) 



(18) 



For all-sky coverage, ±W^/" = ^u'^mm'-i where 5^ is the 
Kronecker delta symbol. Using W , we can simplify (17) to 

±9im{k) = J2[ '^^'^'^ ±WSr'Zee,{k,k')±he,m'{k') (19) 



where 



Zie{k,k') = - / dsfK{s)drp{s\r)j(,{k'r)je{ks). 



(20) 



±htrn{k) may be calculated by direct substitution of the ex- 
pansion of 7± , yielding 

±hem{k) = J dk'k''^Mee{k,k') ±'rem{k') 

where 



(21) 



Mu'{k,k') 



drfjc ir)je' {k'r)je {kr)n{r). 



(22) 



(This equation is more general than necessary at this stage, 
with two indices £ and £', but M^^/ will appear in its general 
form in the shot noise below). It is convenient to define a 
continuous form of the Einstein summation convention in- 
dicating integration over wavenumber: 

A{k,k')B{k',k") = J dk'k'^A{k,k')B{k',k"). (23) 

With this notation, we can write the siiear expansion coef- 
ficients in terms of the present-day, hncar density field coef- 
ficients as follows (from (18), (20), (7), (8) and (6)): 



ji' + 2)! 
{£' - 2)! 



±WSr Zee'{k,k') 



Mu{k',k)^-i^%fi5%^,{k') 
k'^ 



(24) 



where A = 3Qm.H§/{4c^). Note that the sums over £' and 
m' are made exphcitly. 

For all-sky surveys, the average values of ±gim(k) are 
zero; information about cosmological parameters comes from 
their variance (or, more precisely, their covariance). 



2.4 Covariance matrix of ±gim{k) 

Since the present-day hnear power spectrum is a homoge- 
neous field, its covariance matrix is diagonal: 

{5°,^{k)5°r^,(k')) = ^5^(fc - k')5f,,5^^, 



fc2 



(25) 



where 5^ is the Dirac delta function, and Ps{k) is the 
present-day linear overdensity power spectrum. 

The signal part of the covariance matrix is then 



{±gem{k)±gpm.'{k')) = ^^X] 



(f + 2)! 
{!-2)\ 



Note that if coverage is all-sky, the ±W matrices become 
delta functions and the covariance matrix simplifies consid- 
erably to 

{±gim{k)±g*efrn'{k')) = A^'^^^^Zu{k,ki)Zu{k' ,k2) 



Mu{ki,k3)Mu{k2,k4) 



2.5 Shot Noise 



m(k3,kz)'n(.{k4.,k6) „o 



kt 



PKk^)5Xm'-{'i^) 



The shot noise can be computed via standard methods (e.g. 
Peebles 1980), by dividing the volume into cells i containing 
rij = or 1 galaxy. 



±gim{k) = sj^ ^ ni^±ije{ksi) ±2Y("'{ni). 



(28) 



2 ^ 

{ ±gem{k) ±ge>m'{k')) = -'S2{ninjj±i jlj)je{ksi) 

ij 

Ji{ks,) ±2yr*(fti) ±2Yr{nj). (29) 

Assuming the shear field is uncorrelated with the presence or 

absence of a source galaxy, {ninj^'±i j±j) = (7±!7±j). 
Ignoring the small correlations in the smoothed density field, 

(mnj) = (m) {i = j) 

(ni){nj) {i + j). (30) 

The shot noise comes from the term i = j, for which 

2 



:i7±.r) = (i7±.r) + 



(31) 



where the second term will dominate for weak lensing. Thus 
the shot noise may be written 

{ ±gem{k) ±gprn'{k')}sN = ^ J dsfK{s)n{s) 

3t{ks)ji>(k' s) ±Wu^ 



^Mee>{k,k') iWl 



(32) 



using the definition (22). 



3 ESTIMATION OF COSMOLOGICAL 
PARAMETERS 

Various parts of this analysis are dependent on cosmological 

parameters: the matter power spectrum P°(fc); the compo- 
nents of the metric r{z) and /x[r(2)]; the growth rate of 
perturbations g{t). These parameters {{9a}) may be esti- 
mated from the data using likelihood methods. Assuming 
uniform priors for the parameters, the maximum a posteri- 
ori probability for the parameters is given by the maximum 
likelihood solution. For large-scale modes we use a gaussian 
likelihood 



2 In L(g|{6'a}) = constant - det(C) - g • C7 ^ 



• g 



(33) 



where C = 5 + A'" is the covariance matrix, given by signal 
Zff{k. ki)Zf,f(k', k2)Mii{ki,k3)Mii{k2, kA) and noise terms (26) and (32). Note that the average values 
l>'^>)iltU''i, ks) (261°^ ^ ^^^^ °^ ±gem.{k)) is zero, so the information on the 

k^ i5l 5j- (, parameters comes from the dependence of the signal part of 
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C. i.e. we adjust the parameters until the covariance of the 
model matches that of the data. This was the approach of 
Heavens & Taylor (1995) and Ballinger, Heavens & Taylor 
(1995) in analysis of large-scale galaxy data. 

Considering the vector g as the data set has sonic ad- 
vantages over more traditional methods such as the shear 
correlation function or the shear variance. The main one 
is that it is linear in the shear, so the covariance matrix is 
quadratic, and readily calculable. Quadratic estimators such o 
as the shear correlation function have 4th order covarianccs, h 
which can be cumbersome to calculate, even in gaussian, 
linear theory. We do, however, make an assumption that 
the covariance matrix is gaussian, and it requires numerical 
simulations to establish for which ranges of k and (. this is a 
good approximation. 

3.1 Expected errors on cosmological parameters - 
the Fisher matrix 

For a given experimental setup, the Fisher matrix gives the 
best errors to expect, provided that the likelihood surface 
near the peak is adequately approximated by a multivariate 
gaussian. We illustrate the effectiveness of a 3D analysis by 
analysing an all-sky survey with the following details, and 
compare some of the results with a 2D analysis of the same 
data. For smaller surveys, a reasonable approximation is to 
scale the errors by f~i^J^, where fsky is the fraction of sky 
observed. We assume that the average number density is 
n{z) oc 2^exp [-(z/z.)^ '^] (cf (Baugh & Efstathiou 1993)), 
where z, = z„j/1.412 and Zm is the median redshift of the 
survey, which we take to be 1. We assume a source density 
of 30 or 100 per square arcminute, errors on the ellipticity 
of 0.2 or 0.3, and a ACDM model with matter and vacuum 
density parameters of Q,n = 0.3, Qv = 0.7. 

The Fisher matrix is the expectation value of the second 
derivative of the InL with respect to the parameters Oa' 

and the marginal error on parameter 6a is {F^^)aa 
(Tegmark, Taylor & Heavens 1997) . If the means of the data 
are fixed, the Fisher matrix can be calculated from the co- 
variance matrix and its derivatives by (Tegmark, Taylor & 
Heavens 1997) 

Fafi = ^Trace [C-^C,aC-^C,g] . (35) 

For an all-sky survey, this simplifies, because modes with 
different £ and m are uncorrelated {±W^,'^ = Sf^,S^^i), 
so, defining as the covariance matrix for the modes with 
harmonic £ (and different fe), 

Fcf} = \ ^(2£ + l)Trace [{C')-^ C%{C')-' ^g] . (36) 
t 

For illustration, we have considered modes up to £ = 100, 
and 100 k modes spaced equally between k = 0.001 and 
0.1 /iMpc~^. These should be safely in the linear regime. 
The I range can be extended, but at some point the linear 
approximation will break down, and the likelihood expres- 
sion will deviate from a gaussian. The former may be tackled 
by using a nonlinear power spectrum, such as proposed by 



(2D) Fractional Error on Power Spectrum 




Harmonic 1 



Figure 1. Fractional error on the amplitude of the linear mat- 
ter density power spectrum, assuming all other parameters are 
known, from a 2D analysis of the survey, where individual source 
distance information is ignored. The figure shows the improve- 
ment as more harmonics are added up to ^ = 100. 100 ra- 
dial wavenumber modes are considered between k = 0.001 and 
0.1/iMpc-^ Illustrated model is ACDM, with 100 source galax- 
ies per square arcminute, and an ellipticity error of 0.2. Plot is 
cumulative, showing the reduction in error as more modes axe 
included. 



Peacock & Dodds (1996); the latter will need to be assessed 
by computer simulation. 



3.2 Power spectrum amplitude 

Fig. 1 and 2 show the improvement in the error on the frac- 
tional amplitude of the power spectrum as the number of £ 
modes is increased up to ^ = 100. For a 2D analysis, ignor- 
ing the distance information altogether, except for assuming 
n{z) is known, the error is 1.4% for £ < 100. Using the photo- 
metric redshift distance information improves this to 0.9%, 
assuming, perhaps optimistically, that the photometric red- 
shift error is 0.02. 

The improvement depends on the characteristics of the 
survey. Generally speaking, the better the lensing survey 
(higher number density of source galaxies, smaller elliptic- 
ity errors), the better the 3D analysis does in comparison 
with 2D, at least on the scales shown here. The reason is 
simply that the 3D modes are generally noisy, whereas the 
large-scale 2D modes are not. Improving the survey brings 
in more effective 3D modes which have good signal-to-noise. 
The effective number of good radial modes is investigated 
further in the Karhunen-Locve analysis of section 4. 

The Fisher analysis can be extended to calculate the 
errors arising if we estimate band powers from the 3D lens- 
ing data. The Fisher matrix is calculated as before, and its 
inverse used to estimate the errors on the band powers. Fig. 
3 shows what may be achieved realistically for a survey of 
10% of the sky, with a median redshift of 1. Fig. 4 is for an 
optimistic all-sky survey and optimistic errors. 
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1< 1 <100. k<0.100 h/Mpc. fsijy=1.000 








40 60 
Harmonic 1 



100 



Figure 2. As Fig. 1, except that a 3D analysis is used, with 
photometric redshifts having an r.m.s. error of az = 0.02. 




0.03 



0.04 0.06 
k/(h Mpc"') 



0.08 



0.1 



Figure 4. As Fig. 3, except survey is all-sky, intrinsic ellipticity 
dispersion is 0.2. 



10< 1 <100. k<0.100 h/Mpc. fsijy=0.100 




0.02 



0.04 



0.06 
k/(h Mpc"') 



0.08 



0.1 



Figure 3. Band power estimates for a survey of 10% of the sky. 
The source number density assumed is 30 per square arcminute. 
The error bars are marginalised over the other bands. Note that 
the increase in the error at high A; is a result of our truncating the 
analysis at £ = 100. Increasing this would provide more informa- 
tion on the power spectrum at high k. 




Figure 5. The expected error on w for a survey covering 10% of 
the sky to a median redshift of 1, with 100 galaxies per square 
arcminute, and a photometric redshift error of 0.02, and an el- 
lipticity error of 0.3. The plot is cumulative to harmonic £. In 
practice for partial sky coverage the low-£ modes will not be ac- 
cessible, but this will hardly change the results. 



3.3 Vacuum Energy equation of state 

The most exciting prospect of 3D lensing analysis is to mea- 
sure the equation of state parameter, w, of the vacuum en- 
ergy, defined in terms of its pressure and energy density by 
Pv = wpvC^ ■ The equation of state of the vacuum energy 
influences the lensing signal in two ways. Firstly, the growth 
rate of perturbations differs; secondly, the distance-redshift 
relation is changed. I use the growth rate from Abazajian & 
Dodelson (2002) for quintessence models, and the distance- 
redshift relation is obtained from the integral 



where for flat models H(z) = Hq[{1 - r2„)(l -I- z)'^ + 
f2„(l -I- a)3(i+™)]-i/2 jg ^jjg Hubble parameter in terms of 
the present day Hubble parameter i?o, and n„ is the present 
vacuum energy density parameter. 

The sensitivity to w can then be calculated by com- 
puting the Fisher matrix (scalar), equation (36), where the 
derivative of the covariance matrix is obtained by finite dif- 
ferencing. Note that this figure is optimistic in the sense 
that the errors it produces assume that all other parame- 
ters are known, not altogether unrealistic with the success 
of microwave background experiments. On the other hand, 
£ = 40 is a reasonably large scale (~ 5°), and there is more 
information available at higher ^, so 1% accuracy is probably 
achievable. 
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4 KARHUNEN-LOEVE ANALYSIS 



KL modes: L = 10 



KL modes: L = 20 



For all-sky surveys, the coupling of modes is only between 
those of the same i and m, so the data analysis task is 
relatively straightforward, as each set of i and m can be 
analysed separately. For partial sky coverage, the modes 
get mixed, and the covarianco matrix gets very large. Data 
compression techniques such as signal-to-noise eigenmodes, 
Karhunen-Loeve (KL) transformations and the like can be 
extremely valuable in reducing the size of the data set whilst 
having minimal impact on the error bars on the recovered 
parameters. Those have been used extensively in studies of 
large-scale structure and the microwave background (refs). 

The procedure is detailed, for example, in Tegmark, 
Taylor & Heavens (1997), and involves finding linear com- 
binations of the data which are uncorrelated and which arc 
ordered in decreasing order of the information which they 
yield on a parameter. This involves solving a generalised 
eigenvalue problem 



C,ab: 



ACb 





KL modes: L = 30 



10 

Mode number 
KL modes: L = 40 




10 

Uode number 




(38) 



where Oc, is the parameter of interest. The eigenvalues A 
quantify how much information each linear combination 
b • X of the data x provides, and very often this declines 

sharply with mode rmmbor, so the dataset can be com- 
pressed substantially without significantly increasing the er- 
ror bar. Here, we calculate the KL modes appropriate for 
the (natural log of the) amplitude of the matter power spec- 
trum, in which case C,c« = S. The problem then is equivalent 
to finding signal-to-noise eigenmodes, e.g.Vogeley & Szalay 
(1996). Fig. 6 shows the signal-to- noise of the best modes 
for each of the low-order multipoles. We see that there are 
roughly 4 useful modes with S/N > 1. In other words, there 
is some extra information on the power spectrum coming 
from the 3D information, but the added accuracy is limited 
- the analysis had 100 radial wavenumbers. This is broadly 
in agreement with the findings of Hu (1999), who found that 
division of the source population into two, on the basis of 
redshift, improved the error bars on the power spectrum by 
a factor of 2, but that further division yielded little extra 
accuracy. 

It is quite instructive to look at the KL modes them- 
selves. Fig. 7 shows the best radial KL modes corresponding 
with the low-order angular multipoles. We see several inter- 
esting features. Firstly, the best two modes between them 
weight the data reasonably uniformly, apart from a drop- 
off at small distances. This latter behaviour is expected, as 
the lensing is ineffective at low redshift. The near- uniform 
weighting suggests that a 2D survey is not too bad - it corre- 
sponds to equal weighting, and matches quite well the single 
best mode of the KL analysis. For other parameters, this 
may not hold. 

The other interesting and encouraging feature of these 
graphs is that the good modes are long-wavelength, suggest- 
ing that having precise photometric redshifts is not neces- 
sary to get most of the information present in the data. 



5 CONCLUSIONS 

In this paper, I have shown how photometric redshifts may 
be used to perform a full 3D statistical analysis of the shear 



Figure 6. Signal-to-noise eigenmode analysis of the sample sur- 
vey, for a few of tiie low-order angular modes. We see tliat there 
are about 4 useful linear combinations of the radial k modes. 





1000 2000 3000 




1000 2000 




1000 sooo 



Figure 7. KL modes corresponding to the low-order multipoles. 

In each case the best mode is solid, and the dashed, dot-dashed, 

dotted and the dash-dot-dot-dotted lines show the next 4 modes 
in decreasing order of usefulness. Note that a) the best couple of 
modes largely ignore the distance information, as in a 2D survey, 
and b) the best modes are long-wavelength - accurate photometric 
redshifts may not be required. 



field in weak lensing surveys. With photometric redshifts, 
one has an estimate of the lensing shear field at a the (esti- 
mated) positions of the source galaxies in three dimensions, 
and there is no particular reason to throw this away by ig- 
noring the distance information or by dividing the source 
galaxy population into shells. 

Determining the properties of the vacuum energy, which 
makes up ~ 70% of the energy density of the Universe, is 
one the most important current goals in cosmology. These 
properties are not easy to determine accurately using the mi- 
crowave background alone, which has been so successful in 
recent years in pinning down other cosmological parameters 
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(e.g. Spergcl et al. 2003). I show in this paper that 3D lensing 
analysis offers the possibiUty of high-precision measurement 
of the equation of state parameter w (defined such that the 
vacuum has pv = wpvC^), with an accuracy of 1% being a 
realistic possibiUty. With higher angular resolution than in- 
vestigated here, the prospects of getting good constraints on 
w{z) are good, and the possibilities of testing specific pre- 
dictions of w{z) from models of the vacuum energy axe ex- 
cellent. For other parameters, such as the amplitude of the 
matter power spectrum, 3D information reduces the error 
bars by modest factors. A signal-to-noise eigenmode anal- 
ysis suggests that there are a few radial modes which are 
useful for this purpose. This analysis also shows that the 
high signal-to-noise modes have little high-frequency struc- 
ture, so even modestly accurate photometric redshifts can 
be useful. 

In addition to the advantages in reducing the statistical 
error on parameters, photometric redshifts also allow elim- 
ination of a number of possible systematic effects, arising 
from physical processes, which could otherwise eventually 
limit the accuracy of weak lensing studies. For cosmic shear 
measurements in particular, the dominant physical system- 
atic may be the intrinsic alignment of galaxies arising from 
tidal forces during and after formation (Heavens, Refregier 
& Heymans 2000; Croft & Metzlcr 2000; Crittenden et al. 
2001; Catelan, Kamionkowski & Blandford 2001; Mackey, 
White & Kamionkowski 2002; Jing 2002). Heymans & Heav- 
ens (2003) and King &: Schneider (2002) showed how pho- 
tometric redshifts can be used to remove nearby pairs from 
the analysis and essentially remove the intrinsic alignment 
contamination to high accuracy. The advantages in reduced 
statistical and systematic errors make a compelling case for 
obtaining photometric redshifts for cosmic shear surveys. 
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